Magnetic glassy phase in Fei +y Se x Tei_ x single 
crystals 

G Lamura 1 , T Shiroka 2 ' 3 , P Bonfa 4 , S Sanna 5 , F Bernardini 6 , 
R De Renzi 4 , R Viennois 7 ' 8 , E Giannini 7 , A Piriou 7 , N Emery 9 , 
M R Cimberle 1 and M Putti 1 

1 CNR-SPIN and Univcrsita di Gcnova, via Dodccancso 33, 1-16146 Genova, Italy 
2 Laboratorium fur Fcstkorperphysik, ETH-H6nggcrberg, CH-8093 Zurich, 
Switzerland 

3 Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 
4 Dipartimcnto di Fisica and Unita CNISM di Parma, 1-43124 Parma, Italy 
5 Dipartimcnto di Fisica and Unita CNISM di Pavia, 1-27100 Pavia, Italy 
6 Dipartimcnto di Fisica, Universita, di Cagliari, 1-09042 Monserrato (CA), Italy 
7 Departement dc Physique de la Maticre Condcnsce, Univcrsite dc Geneve, 24 Quai 
Ernest- Ansermet, CH-1211 Geneva, Switzerland 

8 Institut C. Gerhardt, UMR 5253, Universite Montpcllicr II, place Eugene Bataillon, 
F-34095 Montpellier, France 

9 ICMPE/GESMAT, UMR 7182 CNRS-UPEC, 2-8 rue Henri Dunant, F-94320 
Thiais, France 

E-mail: tshiroka@phys . ethz . ch 

Abstract. The evolution of the magnetic order in Fei+ySe^Tei-a, crystals as a 
function of Se content was investigated by means of ac/dc magnetometry and muon- 
spin spectroscopy. Experimental results and self-consistent DFT calculations both 
indicate that muons are implanted in vacant iron-excess sites, where they probe a local 
field mainly of dipolar origin, resulting from an antifcrromagnetic (AFM) bicollinear 
arrangement of iron spins. This long-range AFM phase disorders progressively 
with increasing Se content. At the same time all the tested samples manifest a 
marked glassy character that vanishes for high Se contents. The presence of local 
electronic/compositional inhomogeneities most likely favours the growth of clusters 
whose magnetic moment "freezes" at low temperature. This glassy magnetic phase 
justifies both the coherent muon precession seen at short times in the asymmetry 
data, as well as the glassy behaviour evidenced by both dc and ac magnetometry. 
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1. Introduction 

Among the recently discovered Fe-based superconductors, the iron chalcogenide family 
Fei+ySe^Tei-a; is structurally one of the simplest [I], being characterized by layers 
of square lattice of Fe ions, tetrahedrally coordinated with a chalcogen. At room 
temperature its tetragonal cell belongs to the PA/nmm space group, with iron atoms 
occupying also an additional site in the Te plane [2]. The occupation of this additional 
Fe site can be reduced by growing crystals starting from Fe-deficient and Se-rich 
compositions [3]. Tiny variations in both excess iron (y) and Se substitution rate 
(x) strongly affect the structural, electronic and magnetic properties of Fei +y Se x - 
Tei_ x . Upon cooling, the room-temperature tetragonal crystal structure of the parent 
compound (x = 0) becomes either monoclinic (P2/lm), for y < 0.125 [21 Hj, or 
orthorhombic (Cmme), for y > 0.125 [SI E]. This transition to a lower-symmetry 
structure is highly correlated with a magnetic order that involves the Fe moments 
below ~ 70 K. It consists in a bicollinear commensurate antiferromagnetic order, as 
predicted theoretically [7] and proved experimentally by neutron diffraction [H [S]. 
This magnetic order becomes incommensurate when the concentration of excess Fe 
atoms is larger than a critical value (y ~ 0.1076) [H Ej- At the same time, a 
partial substitution of Te with Se progressively suppresses the tetragonal-to-monoclinic 
structural transition, changes the nature of AFM order into an incommensurate short- 
range (SR) one and induces superconductivity at low temperatures. Several neutron 
scattering (SI IS1 QUI UH U21 SI IS1 U31 and zero-field muon spectroscopy (ZF-//SR) 
[T5| [TBI [T7] studies delineate a picture in which the static long-range AFM order evolves 
into a short-range incommensurate one for excess iron values about 0.1 and/or for Se 
substitution rates exceeding 0.1, with the magnetic phase eventually coexisting with 
superconductivity [161 El- This picture is made even more complex by the presence 
of a spin-glass phase, as suggested by dc magnetization data for both x > 0.15 and 
y > 0.01 [HI [HI [191 [20]. Moreover, recent extended X-ray absorption fine-structure 
(EXAFS) measurements as a function of temperature have evidenced local structural 
inhomogeneities in crystallographically homogeneous Fei+ySe^Tei^ samples: Se and 
Te occupy distinct sites with consequently inhomogeneous distribution of the Fe-Se/Te 
bonds [2H [22]. Finally, scanning transmission electron microscopy (STEM) and electron 
energy-loss spectroscopy (EELS) on Fei+^Se^Tei-a; single crystals indicate the presence 
of nanometer-scale phase separation and chemical inhomogeneity [23]. These intrinsic 
inhomogeneities most likely play a major role in determining the fundamental electronic 
structure and the magnetic properties of this compound. 

To date the nature of the magnetic order, crucial for understanding the whole phase 
diagram, is not yet clear and is still under debate. To shed light on the evolution of 
the magnetic ground state as a function of the iron excess and Se doping, we carried 
out a systematic study by means of ac and dc magnetometry measurements as well 
as by zero- and transverse-field /iSR spectroscopy on a series of Fei+ySe^Tei-x single 
crystals. Zero-field //SR represents a unique and powerful local probe for following the 
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Figure 1. (Colour online) Fei+ySe^Tei-a; DC magnetization data for x = 0, 0.07 
(panel a) and 0.21, 0.27 (panel b), respectively. Measurements were carried out in 
an external magnetic field of 200 mT. The vertical arrows indicate the irreversibility 
temperature Tj rr (see text). 



evolution of the magnetic order from long-ranged commensurate, through long-ranged 
incommensurate, down to short-range order, able to distinguish also the presence of 
possible spin-glass/cluster spin-glass phases pU EHJ |26]. By combining this technique 
with complementary magnetization measurements, it was possible to detect the presence 
of a glassy magnetic state induced by Se substitution. After describing the sample 
preparation and characterization in Sec. 12.11 we focus on ac/dc magnetometry and /iSR 
experimental results in the rest of Sec. |2j Successively, we discuss our main findings in 
Sec. El which are then summarized in the conclusions. 

2. Experimental details 

2.1. Sample preparation and characterization 

High-quality, single-phase Fei+ySe^Tei-a; crystals were grown starting from two different 
Fe:(Te, Se) ratios, 1:1 and 0.9:1, corresponding to two different compositions of the 
Te-rich flux. Crystalline samples were grown inside vacuum-sealed quartz crucibles in 
a vertical temperature gradient, by using a modified Bridgman method. The crucibles 
were slowly cooled from the liquid state (T melt = 930-960°C, depending on Se content) at 
cooling rates between 1 and 5°C/h. Further details on sample preparation are reported in 
Ref. [3]. Crystals with actual compositions in the range 0.03 < y < 0.1 and < x < 0.3 
could be obtained. Successively they were thoroughly characterized via structural and 
chemical analysis, as well as by transport and magnetic measurements. The actual 
composition was obtained from single-crystal X-ray diffraction pattern refinement, as 
described in Ref. [3]. Transport properties (resistivity and thermopower) have been 
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reported in Ref. [27]. Selected mm-size crystals of high purity and homogeneity were 
used for the macroscopic (magnetometry) and microscopic (muon-spin spectroscopy) 
magnetic measurements reported in the present work. 

2.2. Ac and dc magnetometry measurements 

The macroscopic magnetic properties of the Fei+^Se^Tei-^ samples were investigated 
via dc and ac magnetometry by means of a superconducting quantum interference 
device (SQUID) and a physical property measurement system (PPMS), respectively 
(both by Quantum Design). In case of dc magnetization we adopted a parallel- field 
configuration, with the sample surface being parallel to the applied magnetic field. 
At low-temperature, low-field magnetization measurements indicate that all the Se- 
substituted samples are superconductors, whose transition temperature and shielding 
fraction increase with Se content, with the latter remaining far from the ideal —1 value 



(see table [Hand Appendix A ). At higher applied fields all samples present interesting, 
albeit complicated, magnetic properties. In figure [1] we show the zero- field cooled (ZFC) 
and the field cooled (FC) low-temperature bulk magnetic response at a static magnetic 
field of 0.2 T. Samples x = and 0.07 display signatures of an AFM order at about 
70 and 40 K, respectively. Superimposed to this magnetic order, there is a hysteretic 
behaviour, as evinced from the bifurcation between the ZFC and FC dc magnetization 
curves. In conventional spin glasses, the temperature below which the ZFC and FC 
curves bifurcate is defined as the irreversibility temperature, T; rr |28j. Interestingly, 
in our case T irr decreases continuously with increasing Se content. To better clarify 
this possible low-temperature magnetic glassy phase, hysteresis cycles were measured 



as well (see Appendix A and figure IA~2) . The cycles show clear "S'-shaped" profiles for 
temperatures below T irr for almost all the samples. Measurements of ac susceptibility 
were carried out on all the samples under test at three different frequencies (0.1, 1, and 10 
kHz), with an ac excitation field H ac = 1 mT. No clear evidence of frequency dependence 
was found, except for the x = 0.21 sample. Figure [2] shows the low-temperature ac 
susceptibility data Xac(T). At the lowest frequency, the temperature where Xac shows 
a maximum is generally defined as the freezing temperature, Tj [28]. We remark that 
this is similar to the T ur of the dc measurements. As the frequency of ac measurements 
increases, we observe a continuous increase in Tf. 

2.3. Muon-spin spectroscopy 

Muon-spin spectroscopy (/iSR) experiments were carried out at the GPS spectrometer 
(7rM3 beam line) of the Paul Scherrer Institut, Villigen, Switzerland. A set of five single 
crystals was measured both in transverse- (TF) and in zero-field (ZF) configurations. 
The technical details regarding the sample holder and the measurements are reported 
in Appendix B A schematic drawing of the spectrometer, including the sample 
orientation, the beam direction and the detector positions is shown in figure [3j TF- 
/iSR measurements were carried out by applying a (small) magnetic field perpendicular 
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Figure 2. (Colour online) AC susceptibility of Fei+ySe^Tei-a for y ~ 0.049 and 
x = 0.21. Inset: frequency dependence of T/, the temperature of % ac maximum, as 
measured at v = 0.1, 1, and 10 kHz (see Sec.[3]for details). 
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Figure 3. (Colour online) Geometry of the ZF-^SR experiment showing the detector 
positions, the sample orientation, as well as the directions of muon momentum and 
the initial muon-spin polarization (normal and rotated). 



to the muon momentum p^. This configuration provides a straightforward way for 
estimating the paramagnetic volume fraction from the loss of spin polarization in the 
complementary magnetic part (see below). ZF-/iSR measurements were carried out 
using both a rotated and a nonrotated muon-spin polarization. In both configurations 
the samples were oriented so as to have the direction of the incident muon beam 
orthogonal to the sample surface and, therefore, parallel to the c axis of Fei+ySe^Tei-a; 
crystals (see figure [3]). 

2.3.1. Transverse- field fiSR All the selected crystals were investigated in a transverse 
field of 5 niT at two different temperatures: well below and above the respective Neel 
temperatures, T^. The latter were determined by dc magnetization measurements and 
confirmed also by ZF-yuSR (see next paragraph). In an applied transverse field, the 
coherent precession of muon spins in a nonmagnetic environment gives rise to a constant 
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Figure 4. (Colour online) Low-temperature ZF-/1SR asymmetry as detected in the 
x = 0.07 case by the F-B and U-D detectors in a configuration where the muon spin 
and momentum are collinear (and perpendicular to the iron planes). 

oscillation amplitude. On the other hand, muons stopping in magnetically ordered parts 
of the sample sense a wider distribution of fields, with a quick loss of coherence and 
consequent reduction (or even "disappearance" ) of the relevant asymmetry signal. Due 
to the random character of muon implantation, the amplitudes of the different types of 
precession signals are proportional to the volume fraction of the respective environments. 
Hence, by evaluating the ratio of the transverse asymmetries measured at T <C Tn and 
at T Tn we could determine the paramagnetic volume fraction, Vp ara /Vtot, as reported 
in table [TJ 



Table 1. Magnetic properties of various Fei+ySe^Tei^a; single crystals as determined 
from fiSK and ac/dc magnctomctry measurements. 
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2.3.2. Zero- field jj,SR ZF measurements were carried out using either a rotated or 
a nonrotated muon-spin configuration, depending on the system availability. Typical 
asymmetry spectra for the x = 0.07 case, collected at two selected temperatures in the 
forward-backward (F-B) detectors with a nonrotated muon spin, are shown in figure HJ 
One can distinguish two different temperature regimes: in the low-temperature regime 
(T < Tn = 43.6 K) clear oscillations with a rather high damping rate are observed, 
demonstrating that the sample has made a transition to a magnetically ordered state 
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characterized by long-range order. Here the oscillating signal is superimposed to a slowly 
decaying exponential relaxation, observable only at longer times (not shown). Above 
In, instead, the amplitude of the oscillating signal goes abruptly to zero, leaving only 
two slow but distinct exponential relaxations. To follow the above evolution across the 
whole temperature range, all data sets were fitted to the function: 

A(t) =a T -f (B„, t) ■ exp(-4t 2 / 2 ) + a L • exp(-At) + a bg . (1) 

Here the transverse component, characterized by an amplitude a-r and a damping rate 
o"t, is an oscillating function f(B^,t) = J (7 M -B M t), with J the zeroth-order Bessel 
function, the magnitude of the internal magnetic field, and 7^ = 2ti x 135.5 MHz/T 
the muon's gyromagnetic ratio. The longitudinal component was fitted to a relaxing 
Lorentzian term, with aL and A ~ 0.1 /is -1 , the initial amplitude and relaxation rate, 
respectively, while db g is a temperature-independent relaxing background term that 
could be ascribed to a non negligible volume fraction that remains paramagnetic even 
at low temperatures (as inferred also from TF-/xSR data — see previous section and 
table d]). It is worth noting that the highly damped oscillation found at short times 
below Tn in the x = sample could be fitted satisfactorily using either a cosine term or a 
Bessel function. Despite the long data-acquisition times, the collected statistics was still 
not high enough to clearly distinguish between these two cases. Consequently, although 
we can safely state that the system shows a long-range magnetic order, not much can 
be said about the commensurability of the order [T7]. By taking into account the whole 
/iSR data set and the theoretical predictions on the muon implantation sites (vide infra), 
for the x = case (pure FeTe) we privilege the fit results where the oscillating term 
was fitted with a zeroth-order Bessel function. This represents a coherent choice, in 
agreement with the results for the x = 0.07 and x = 0.21 samples, where the best 
fits were obtained with a Jo term, too. As for the sample with x = 0.27, here the 
damping was too strong and no oscillations could be detected at all [f(t) — 1]. In 
figure |5] we summarize the temperature dependence of the amplitudes of the transverse 
component ay. Since a transverse component in ZF data can arise only in presence 
of a magnetically ordered phase, by fitting the temperature-dependent amplitude czt 
with a complementary erf function, 1 — erf[(T — Tn)/(v / 2ATn)], we could determine 
the Neel temperature, Tn, and the transition width, ATn, of the ordered phase. The 
resulting values for each sample are reported in table [U and in figure EJ We note that 
T N decreases progressively as the Se content increases. Figure [7] shows the temperature 
dependence of the internal magnetic field B^, as probed by the implanted muons. The 
best fits at low temperature yield a mean-field value comprised between 100 and 200 
mT for samples with x — 0, 0.07, and 0.21 (right scale in figure [7]). In case of x = 0.27 
the oscillations disappear and, in absence of £>„, the width of the Gaussian relaxation, 
(7/7^ = (B 2 — B^ yi 2 [291 EQ], can only provide an upper limit to the distribution of 
the internal fields at the muon site. In our case, at the lowest measured temperatures, 
these field values lie in the range 50 to 150 mT (left scale in figure [7]). 
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Figure 5. (Colour online) Temperature dependence of ar, the amplitude of ZF-/iSR 
transverse component, as measured in different Fci+^Se^Tei-^ single crystals. In 
each case the value of ot is proportional to the magnetically ordered volume fraction 
(maximal at T = 0). All the data has been collected by the U-D detectors, except 
for the x = 0.07 case, where no spin rotator could be used and the data refer to the 
forward-backward (F-B) detector pair. Solid lines represent numerical fits with an erf 
function (see text for details). 




x(Se) 



Figure 6. (Colour online) Tn versus Se content for all the tested samples, 
as determined by fitting the temperature dependence of the ZF-/xSR transverse 
component with an erf function. The line is a guide for the eyes. 

2.3.3. Direction of the internal magnetic field In order to study the direction of the 
internal magnetic field in the differently doped samples we take advantage of the 
ZF-/iSR as a non-perturbing technique (no applied external field). To simplify the 
interpretation of the results we consider here the case when the initial muon-spin 
polarization is antiparallel to the beam direction and systematically compare the data 
collected by the F-B with those of the U-D detectors (see figure [3]). We recall that 
each investigated sample consists of a mosaic of single crystals from the same batch, 
positioned so that their ab surfaces coincide with the xy plane (i.e. orthogonal to the 
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Figure 7. (Colour online) Temperature dependence of the internal magnetic field 
(□) and of the transverse relaxation Gaussian width, (• ), for the investigated 

Fei+ySe^Tei-a; single crystals. 



muon spin), as shown schematically in figure |3j Therefore, the initial spin polarization is 
always parallel to the crystalline c axis, while the a and b axes of the mosaic crystals are 
randomly oriented in the xy plane. As an illustrative example, figure H] compares the 
muon-spin asymmetry measured in the up-down vs. forward-backward detector pairs 
for the x = 0.07 case at T = 5 K. Similar results were obtained also for other samples 
and/or different temperatures (not shown). We note that the asymmetry measured 
by the U-D detector pair (i.e. along the x direction) is practically zero||| At the same 
time the signal detected by the F-B pair (i.e. along the z direction) not only accounts 
for the total measured asymmetry, but its value A to t(t) changes with time, reflecting 
a magnetically ordered state. These features indicates that the internal magnetic field 
B^, as probed by the implanted muons, can be exclusively either in the ah plane or 
orthogonal to it (i.e. not in some arbitrary out-of-plane direction). Now, if the internal 
magnetic field _B M were parallel to the crystal c-axis, we would expect the asymmetry 
measured by the U-D and F-B detector pairs to be equal to zero and to a constant value, 
respectively. Since this is not the case, we conclude that the direction of the internal 
field lies necessarily in the laboratory xy plane, i.e. in the sample's ah plane. 

\ Supposing Bfj, to be in the sample's a&-plane, a small c-axis misalignment (6 C < 10°) could account 
for a tiny nonzero asymmetry in the U-D detectors: oud oc ofb • tan# c ~ 0.02 independent of 
time / temperature . 
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3. Discussion 

We discuss now the experimental results presented so far with the aim of obtaining a 
unified picture for the entire Fei+^Se^Tei-a; series. 

DC magnetization: as a general feature, the investigated samples display a complex 
magnetic behaviour. First, all the samples with x > are superconductors, whose T c 
and shielding fraction scale with x. Since the normalized SC shielding is always far 
away from the ideal value —1, even for the sample with the highest Se content, this 
suggests that superconductivity has not a bulk character. Second, x(T) curves taken 
at 200 mT (large enough to destroy a weak superconductivity) exhibit a pronounced 
hysteresis, with irreversibility temperatures T irr that decrease both with increasing Se 
content (instead of being independent of it, as previously reported, e.g., in Ref. [H]) and 
with the applied field (not shown). In particular, the field dependence of T irr is generally 
ascribed to the presence of a spin glass phase [28]. Such a hysteretic behaviour most 
likely originates from the competition between the antiferromagnetic mtra-cluster and 
the ferromagnetic mter-cluster interactions. It suggests the lack of a true long-range 
magnetic order (STJ [32j [33j EH [35] , as confirmed also by our ZF-//SR experimental results 
(vide infra). Further insight on the glassy magnetic behaviour is obtained from the 
isothermal magnetization curves taken below and above the respective antiferromagnetic 
transitions (see figure lA2l) . The low temperature data show two noteworthy features: 
i) for x(Se) > the first magnetization curve is compatible with the presence of a small 
superconducting phase, ii) The hysteresis cycles consist of a complicated superposition 
of weak superconductivity (for x > 0) and a glassy magnetic state, that give rise 
to characteristic " 1 S'-shaped" m(H) profiles (281 EH1 E7]. This feature becomes more 
prominent at a higher Se content, whereas it seems hidden in the x = 0.07 sample. 
The close lying superconducting and irreversibility temperatures make it very difficult 
to distinguish between these two coexisting effects. 

AC susceptibility: only the ac susceptibility data of the x = 0.21 sample display a clear 
frequency-dependent maximum at low temperature. Figure [2] shows the temperature 
dependence of Xac, the ac magnetic susceptibility. The maximum of Xac, observed at 
freezing temperatures Tf = 14-15 K for frequencies v = 0.1-10 kHz, represents the 
point where the system's relaxation time r equals the characteristic time scale of the 
ac measurement, 1 jv. The ac magnetic response becomes frequency- independent only 
above that maximum, a behaviour peculiar of spin glasses [28] . The Tf value determined 
from the Xac maximum is in very good agreement with the T irr value extracted from the 
dc-magnetization FC-ZFC irreversibility curves at 200 mT (see figure []]). We notice 
also that the peak position, Tf, shifts to higher temperatures with increasing frequency, 
yet another characteristic of spin glasses or superparamagnetic states. In fact, in these 
systems the relative freezing temperature variation, ATf/Tf, scales proportionally with 
the logarithmic frequency variation of the excitation frequency A log(i') [381 Ell HQ]: 

ATf/Tf = 5 ■ A \og(u) (2) 



Magnetic glassy phase in Fe\ +y Se x Te\- X single crystals 



11 



The scaling parameter, 5, usually lies between 0.005 and 0.06 in case of spin glasses, 
while it exceeds 0.1 in superparamagnets with noninteracting spin clusters. As shown 
in the inset of figure El the fit of experimental data provides 8 = 0.029(1), a value close 
to those found in other spin- glass/cluster- glass systems [38] EJ H21 [28]. This rules out 
any possible superparamagnetic state and represents strong evidence about spin-glass 
behaviour in the Fei+ySe^Tei-a; family. 

Muon spectroscopy: the low-temperature, short-time asymmetry data for < x < 0.2 
show a coherent muon-spin precession successfully fitted by a Jo term. This feature is 
generally ascribed to the onset of a long-range AFM order, although somehow weakened 
by the incommensurate character of the magnetic structure [13] or, more simply, by the 
presence of magnetic clusters (HI 121] . In principle, our findings are in agreement with 
previous neutron diffraction [6] [TT] [4] [5] and /xSR measurement [16] [17] data, both sug- 
gesting a long-range AFM magnetic order. However, the more marked manifestations of 
a disordered magnetism in our case, lead us to associate the Bessel term of the fits with 
a generic medium-to-short range magnetism, as is the case of glassy magnetic phases. 
Indeed, the presence of large magnetic clusters has been invoked to explain the damped 
oscillations observed in certain conventional spin-glass systems [261 [25] . This picture 
most likely applies whenever sufficiently big magnetic clusters can cause both coherent 
muon precessions, and give rise to a glassy magnetic hysteresis in the dc magnetization 
data. As the Se content is further increased beyond x ~ 0.2, there is a progressive loss 
of muon-spin polarization, reflected in an enhanced damping of the oscillating signal, 
which evidences a transition to a pure short-range magnetic order [13]. Further hints 
towards the presence of a glassy magnetic phase can be inferred from the temperature 
behaviour of the longitudinal relaxation Ai (see figure El). Here the broad peak that 
develops below Tn can be explained in terms of a progressive freezing of the magnetic 
correlations. In general, magnetic correlations with a glassy character were found in all 
the samples under test. Indeed, the typical features expected in a classical spin-glass 
phase (such as, e.g., the bifurcation in dc magnetization) are present in all of them, 
but they are most pronounced in the x = 0.21 case, which fulfills at best the criteria 
defining the spin-glass state (see also Figs. 121 and |A2]) . In this sense, by considering 
only dc magnetization data, our results fully agree with the phase diagram reported in 
Ref. [H]. However, by considering also the additional data of ac susceptibility, muon 
spectroscopy and dc isothermal magnetization, we obtain a scenario where Fe^Se^- 
Tei_ x seems to display a generalized magnetic glassy phase, characterized by hysteretic 
magnetic correlations from x = up to x = 0.27. As detailed in our final remarks, the 
difficulty in unambiguously defining the spin-glass state in Fei+ySe^Tei-a; is most likely 
due to the presence of local structural, chemical and electronic inhomogeneities, intrin- 
sic to this class of materials [21] [22] [23]. As described in Sec. I2.3.3[ our experimental 
data suggest that the internal magnetic field, as sensed by the precessing muons, lies in 
a plane parallel to the iron layers. To better appreciate the implications of this finding 
we tried to determine the muon-stopping site(s) by investigating an ideal stoichiometric 
FeTe crystal with the Density Functional Theory (DFT) formalism in the generalized 
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Figure 8. (Colour online) Temperature dependence of the muon-spin longitudinal 
relaxation rate A for the investigated Fci+ySc^Tei-z single crystals. The light blue 
areas indicate the magnetic phase. 

gradient approximation (GGA) |45]. We adopted the Full Potential Linearized Aug- 
mented Plane Waves (FP-LAPW) method, as implemented in the Elk code [2B]. In 
particular, we considered a bicollinear magnetic order on a 2a x a x c cell with PA/nmm 
symmetry [6] and a -R^in x niax(|k|) = 8 is chosen for the expansion of the wave func- 
tions in the interstitial region (-R^in * s the smallest muffin-tin radius in the unit cell and 
max(|k|) the largest wave number of the basis set). The resulting magnetic moment per 
iron atom was 2.24 /xb, in reasonable agreement with 2.54 ub, the experimental value 
obtained via neutron powder diffraction [6]. To determine the muon-stopping sites we 
considered the minima of the sign-reversed electrostatic potential, as obtained from a 
self-consistent DFT calculation. The positively charged muons will thermalize in the 
interstitial positions corresponding to the potential minima. Two inequivalent minima 
were found, whose positions are shown in figure M, with the respective coordinates re- 
ported in table [2j The absolute minimum (site A) replicates the one determined in 
Ref. [47] for FeSe, while the secondary minimum (site B) coincides with the centre of 
the tetrahedron formed by Te atoms. The latter site turned out to be unstable against 
the zero-point motion of the muon. This was accounted for by modelling the potential 
minimum with an anisotropic harmonic well (AHW) V(r) = ^m(u 2 x 2 +u 2 y 2 +u 2 z 2 ) and 
evaluating the zero-point energy E of the implanted muon. The resulting isosurfaces 
for V A ' B (r) = V^i'jf + Eq' b are presented in figure fTUl While A-site minima remain on 
disconnected isosurfaces, those for B sites form an interconnected network that rules out 
the presence of muon bound states. Therefore, also those muons initially implanted in 
sites of type B will migrate towards the nearest A site. Hence, sites of type A represent 
the only possible muon implantation site in the FeTe crystal. We recall that sites A 
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Figure 9. Possible muon stopping sites within the tetragonal unit cell (high- 
temperature paramagnetic phase), as obtained from Coulomb potential minima. The 
positions of Fe and Te atoms were taken from Ref. [6]. The arrows represent Fe-ion 
magnetic moments in the ordered phase. Upon considering the fi + zero-point motion, 
only sites of type A are found to be stable. 




Figure 10. Isopotcntial surfaces V(r) for sites A and B calculated at the respective 
V A ' B (r) = V 7 2^+E A ' B (see text). The disjointed inner surfaces (red) refer to A-sites, 
while the interconnected sheets (yellow) to the B-sites. 

Table 2. Properties of the muon-stopping sites. The reported positions correspond 
to the minima of the potential V(r). The respective Wyckoff sites are A 2c and B 2b. 
Eq represents the zero-point energy (see text) and WFC is the wave function centre. 



Site 


Position 




Zfl+ WFC 


V-V A (cV) 


E (cV) 


A 


(V4, V4, 


0.245) 


0.27 





0.36 


B 


(3/4, 1/4, 


0.500) 


0.50 


0.9 


0.39 



almost coincide with the positions of Fe(II) but, since their iron occupancy is lower than 
10% for all the considered y compositions (see also table [Hand Ref. [3]), they are widely 
available for hosting the implanted muons. Note that the electrostatic potential has 
an anharmonic component along the c direction. Therefore, the equilibrium position 
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of the muons may not correspond to the potential minimum. Indeed, considering the 
anharmonicity, we find that the correct value for the site A is ( 1 /4, 1 /4, 0.27), as obtained 
from the centre of the ground-state wavefunction. Taking this into account, we can eval- 
uate the local field probed by the implanted muons, B^, which consists mostly of two 
contributions, Ba and B c , the so-called dipolar and Fermi-contact fields, respectively 
|43j . If we assume that the behaviour of Fei+^Se^Tei-a; series does not differ much from 
that of pnictides, then the Fermi contact contribution is expected to be negligibly small 
|48j . The dominant Bd term at the muon-stopping site was calculated by means of the 
probability density obtained from the AHW ground-state wavefuction. By considering 
the ordered configuration of Fe(I) magnetic moments (each carrying 2.54 /ib, as from 
neutron powder diffraction on Fei +?/ Te polycrystalline samples) j6], the expected dipolar 
field at the muon implantation sites is ~ 238 mT. The latter result is in good agreement 
with 197 ± 11 mT, the value obtained from ZF-//SR data. The residual discrepancy can 
be attributed to a possible small contribution arising from the Fermi-contact term, that 
was neglected in a first instance. We recall here that the excess of iron (present in all 
samples) is strongly magnetic [19] and may play an important role in broadening the 
local field distribution. Moreover, as shown in figure U\ the low-temperature field value 
at a muon site, decreases slightly with increasing Se content. Provided the hyperfine 
coupling tensor does not depend (or depends only weakly) on Se content, this would 
imply that the magnetic moment of the iron ions does the same, in apparent contrast 
with neutron powder diffraction results [6]. 

Finally, we recall that EXAFSJUJ [22], STEM and EELS [23] results, previously men- 
tioned in the introduction, bring strong and clear evidence about local structural inho- 
mogeneities, nanometer-scale phase separation and chemical inhomogeneities on crystal- 
lographically homogeneous Fei+ySe^Tei.^samples. These important findings let us unify 
our experimental results into a coherent picture. Since X-ray diffraction measurements 
evidence that our samples are in a single phase [3] , the above mentioned nanoscale phase 
separation scenario would imply: i) the formation of magnetic clusters, hence justifying 
the glassy behaviour evidenced by both ac/dc magnetometry and muon spectroscopy 
data. This same glassy behaviour could potentially attenuate, or even mask, the true 
variation of the Fe magnetic moment with Se/Te substitution, ii) The presence of a non 
negligible volume fraction that stays in the paramagnetic state also at low temperatures. 
In our case, this fraction represents ca. 20-30% of the sample volume and is mostly Se- 
independent. Hi) The absence of any possible nanoscopic-like M-SC coexistence, rather 
common instead in the 1111 pnictide family [50j ED, [52]. The latter can be ruled out 
because none of the investigated SC samples shows bulk superconductivity. In partic- 
ular, this seems to be a general feature of the iron chalcogenide family that confirms a 
posteriori previous muon spectroscopy and neutron diffraction results on Fei+ySe^Tei-x 
single crystals [TT] IT6]. 
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4. Conclusion 

The low-temperature antiferromagnetic order in an extended series of Fei+ySe^Tei-a; 
single crystals was studied by means of ac/dc magnetometry and by muon-spin 
spectroscopy. Muons are most likely implanted close to a Fe(II) site, where they sense a 
local magnetic field ~ 200 mT, lying in the ah plane. At low Se content /iSR detects 
a medium/long range AF order that progressively weakens with Se substitution and, 
for x > 0.2, turns into a short-range order. Both ac and dc magnetometry data on the 
same samples show clear evidence about a magnetic glassy phase, which becomes more 
prominent for x = 0.21, but vanishes at higher Se content. The magnetometry results, 
combined with the unusually highly-damped ZF-/iSR oscillations, observed in the same 
temperature range where the glassy behaviour occurs, suggest that the system is most 
likely in a magnetic cluster-glass state, probably induced by the nanometer-scale phase 
separation and the chemical inhomogeneity. The last feature seems to be an intrinsic 
property of the Fei+^Se^Tei-a; family, as confirmed also by recent extensive studies on 
polycrystalline FeSe .5Te .5 [53] . 
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Appendix 

Appendix A. Supplementary information on DC magnetization 

We present here supplementary material on dc magnetization measurements. It consists 
in: a) the low-field temperature dependence of the dc magnetization, reported in 
figure lATj and b) the hysteresis cycles measured below and above T N , as determined by 
ZF-/iSR spectroscopy, reported in figure IA21 

Appendix B. The /iSR sample holder and measurement details 

The /iSR measurement of single crystals with a typical size of 3 x 4 mm 2 and 0.1- 
0.3 mm thickness was made possible by using a standard "mosaic technique". Several 
crystals from the same batch were sticked close together using an aluminated mylar 
tape. The resulting mosaic was then centred on a standard aluminium fork sample 
holder suited for muon fly-past measurements. For all the investigated samples the 
effective sample area exposed to the beam was never below 25 mm 2 . Taking into account 
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Figure Al. (Colour online) Shielding volume fraction for all the Se-substitutcd 
samples obtained from ZFC curve as measured by dc SQUID. The data were taken 
in an applied magnetic field of 0.5 mT (except for the x = 0.07 sample, which was 
measured at 10 mT). The vertical arrows indicate the onset of superconductivity. 





Figure A2. (Colour online) Hysteresis cycles at 5 and 100 K for all the Fei+ySe^Tci-a; 
samples under study. Insets show an expanded view of the magnetic field dependence 
of the first magnetization curve at 5 K. 
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the relatively large sample thickness (> 0.25 mm) and the average sample density, a 
0.6 mm thick kapton mask was sufficient to ensure that the majority of the incoming 
muons were implanted into the sample volume. The further use of veto counting (which 
accounts for only those muons effectively stopped in the sample), contributed also to a 
significantly improved signal-to-noise ratio. In this configuration a total asymmetry in 
the range 0.15-0.22 could be obtained, still low if compared with the normal asymmetry 
value expected at the GPS instrument (~ 0.28). The discrepancy could be ascribed to 
muonium (or muonated radical) formation in the kapton mask with a relatively high 
hyperfine coupling constant. Since the ensuing muon-spin depolarization reduces only 
the total signal asymmetry, but does not otherwise affect the measured signal, we could 
safely ignore the fraction of muons being stopped in the kapton frame. 
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